Correlation function for the one-dimensional extended Hubbard model at quarter filling 
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We examine the density-density correlation function in the Tomonaga-Luttinger hquid state for the one- 
dimensional extended Hubbard model with the on-site Coulomb repulsion U and the intersite repulsion V at 
quarter filling. By taking into account the effect of the marginally irrelevant umklapp scattering operator by 
utilizing the renormalization-group technique based on the bosonization method, we obtain the generalized 
analytical form of the correlation function. We show that, in the proximity to the gapped charge-ordered phase, 
the correlation function exhibits anomalous crossover between the pure power-law behavior and the power-law 
behavior with logarithmic corrections, depending on the length scale. Such a crossover is also confirmed by the 
highly-accurate numerical density-matrix renormalization group method. 

PACS numbers: 71.10.Fd, 71.10.Hf, 71.10.Pm 



I. INTRODUCTION 



One-dimensional (ID) electron and spin systems have been 
attracted much attention since they often exhibit nontriv- 
ial quasi-long-range-ordered behavior due to the large low- 
dimensional quantum fluctuation effects.-i^ The critical be- 
havior in the ID systems, which is called the Tomonaga- 
Luttinger liquid (TLL) state, has a long history of research 
and the low-lying modes are known to be described by collec- 
tive gapless excitations and physical quantities show power- 
law behavior in the temperature and/or distance dependences. 
It has also been recognized that, in the systems with spin- 
rotational symmetry, logarithmic singularities appear in the 
magnetic -field-dependent corrections to the magnetization^ 
and the spin susceptibility,^ and in the temperature-dependent 
corrections to the spin susceptibility,^ specific heatfi and nu- 
clear magnetic resonance,- etc. Motivated by such develop- 
ments in the field, a number of numerical studies has been suc- 
cessfully performed to examine the logarithmic corrections in 
spin-chain systemS] ^'^°-'^-'^'^^-'"^-'^'^^ The spin-spin correlation 
function of the ID 5* = i Heisenberg model has been exten- 
sively studied as a most fundamental theoretical model to ex- 
amine the presence of logarithmic corrections ! '^•'^i'^'^" At this 
moment, the correlation amplitudes in the asymptotic form 
of the correlation function can be exactly obtained . ^''^^'^-^ In 
contrast, only few efforts have been devoted to those of the 
Hubbard model due in part to the difficulty of analyzing.^"* 
Therefore, the situation is much less satisfactory as far as the 
logarithmic coiTections in the Hubbard model are concerned. 



quarter-filling is given by 



H 



^U^n^^^nj^l + V^njUj+i, (1) 



where ,, is the annihilation operator on the jth site with spin 



s(=t, j), and the density operators are nj 



,T + and 



The hopping energy between 
the nearest-neighbor sites is represented by t. It is known that, 
at zero temperature, the gapped charge-ordered (CO) state 
emerges in the large repulsive U and V region, where the 
phase boundary is determined numerically by using the exact 
diagonalization2^ and the highly-accurate density-matrix 
renormalization group (DMRG) method.— The mechanism of 
this quantum phase transition has also been addressed by the 
bosonization technique and the renormalization-group (RG) 
methodia22i^ In this paper, we perform the detailed analy- 
sis on the equal-time correlation function N{x) = {ujUj^x) 
in the TLL phase. The exponent of the correlation functions 
is characterized by so-called the TLL parameter Kp,~* Espe- 
cially by focusing on the correlation function near the bound- 
ary to the CO insulating state, we show that it exhibits the 
nontrivial crossover, depending on the length scale, from the 
power-law behavior with logarithmic correction for short dis- 
tance, to the pure power-law behavior for large distance. 



In the present paper, we focus on the logarithmic correc- 
tions in the equal-time density-density correlation function of 
the quarter-filled Hubbard model including the Coulomb re- 
pulsion between electrons on site U and the nearest-neighbor 
sites V. So far this model has been analyzed as a minimal 
model to describe physical phenomena in organic solids,-^ 
The Hamiltonian of the ID extended Hubbard model at 



The present paper is organized as follows. In Sec. II, the 
analytical form of the correlation function is obtained by uti- 
lizing the RG technique based on the bosonization method. 
In Sec. Ill, the analytical results are confirmed by using the 
highly-accurate DMRG method. The summary is given in 
Sec. IV. Detailed derivation of the analytical form of the cor- 
relation function is given in the Appendix. 
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II. BOSONIZATION APPROACH 

In this section, we derive the generahzed analytical form 
of the correlation function in the TLL state. We analyze the 
?7 — > oo limit case and the finite U case separately, since 
the picture of the U ^ oo can become transparent with the 
analogy of the spin-chain system which properties are well 
understood. 



A. The U ^ oo limit 

In the U ^ oo limit, since the double occupancy of elec- 
trons is excluded, the extended Hubbard model [Eq. ([T])] re- 
duces to the spinless half-filled model: 



H.c. 



(2) 



where 



It is well known that this model 



_ 1 

"3 ^ "'3"'3 2- 

can be mapped onto the XXZ spin-chain model by using the 
Jordan- Wigner transformation^'^! ^jj^j ^jje physical properties 
have been extensively studied with both the exact treatment 
based on the Bethe ansatz and numerical approaches. Here 
we examine the analytical form of the correlation function by 
using the exact results obtained in the context of spin-chain 
problems. 

The density operator is expressed in terms of the bosonic 
field operator ip a a^-'^'-^' 



It is well-known that the TLL phase is realized for small 

V < Vc{~ 2t), while the gapped CO state appears for 

V > Vc- In the TLL phase, the parameter 7] is known ex- 
actly from the Bethe ansatz as 77 = 1 — tt^^ cos^^{V/2t), 
and the TLL parameter Kp = l/(47y) varies within the range 
of 1/4 < Kp < 1/2 for < < Vc. The TLL param- 
eter Kp approaches to an universal value Kp = 1/4 when 

V ^ Vc- The mechanism of this quantum phase transition 
has also been addressed by the bosonization technique and 
the RG metho d^'^^'-^" and it has been clarified that the 1/4- 
filled umklapp scattering has a crucial role in making the TLL 
state into the gapped CO state. It has also been shown that 
the universality class of this TLL-to-CO phase transition is 
in the Kosterlitz-Thouless transition, where the umklapp scat- 
tering is irrelevant in the TLL phase while it becomes rele- 
vant in the CO phase. It is worthwhile to note that, on the 
phase boundary between the TLL and CO states, the umklapp 
scattering term becomes marginally irrelevant, where it shows 
very slow scaling-parameter dependence and can give rise to 
anomalous corrections to physical quantities. We analyze this 
effect on the density-density correlation function by using the 
RG method developed in Ref. [l^ The resultant form of the 
correlation function for [/ = 00 in the TLL state is given by 



N{x) 



-A 



cosAkpx [1 — {a/x) 



[1 + ia/x) 



2-8K, 



-,3/2 



(4) 



p{x) = 



a 



1 dip 
27r77 dx 



i-iy 



(3) 



where a is the lattice constant and we will set a = 1 in the 
following. The parameter rj can be related to the TLL param- 
eter by Kp = 1/(477). The nonuniversal parameter c will be 
shown later The model Hamiltonian can be expressed in 
terms of the bosonic field ip. 



where kp is the Fermi wave vector, kF{~ 7r/4). The deriva- 
tion of this form is given in the Appendix. The parameter a 
is the short-distance cutoff which appears in the RG method. 
The coefficient A2 is given by A2 = c^/(27r^) by using the 
parameter c in Eq. (O. In Refs.|2l] and|23 the exact form of 
the parameter c has been proposed, where it is controlled by 
the TLL parameter, i.e., c — c{r]), and its exact form is given 
by 



J 



0(77) 



r(77/(2 - 277)) 

2V^r(l/(2-27?)) 



- l/2r/ 










cxp 


2 





sinh[(277 - l)t\ 
t ysinh(77t) cosli[(l — ?/)<] 



r 



277-1 



-2t 



(5) 



The numerical values of A2 are shown in Table. I. For 
([/, V) = (cx), 0), the quantity c{rj) becomes c{rj = 1/2) = 1. 
We note here that the quantity a in Eq. (|4]i is the only unknown 
parameter, which is to be determined numerically. 

Here we find that equation ^ has two different asymp- 
totics; (i) In the short-range region, the power-law behavior 
with logarithmic correction is obtained, while (ii) the logarith- 
mic correction disappears in the long-range region. By not- 



ing (x/a)- 



-{8K, 

which characterizes crossover between these two 



exp[-(8Xp - 2) ln{x/a)], the length 



regions is given by 

Across = ac:xjp[l/{8Kp - 2)]. (6) 

We note that, for V = Vc{= 2t), the logarithmic correc- 
tion appears in the whole length scale. By noting A2 
(2 - V/t)-^/^/{V2TT) for V -> Vc, the expHcit form of the 
correlation function ^ atV = 2tis given by 

1 1 cos4:kFX 



N{x) 



47r2a 



(27r)3/2 : 
for ([/, V) 



■\n^/\x/c 
(00, 2t). 



(7) 
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TABLE I: Amplitude A2 = (r;)/(27r^) as a functions of V/t for U ^ 00. 



V/t 




1//-/- 
V /t 


A„ 
A2 


V/t 




1//-/- 
V /t 


A2 


V /z 


A2 


0.0 


0.0506606 


0.5 


0.0741643 


1.0 


0.107134 


1.5 


0.164769 


1.95 


0.404613 


0.1 


0.0548623 


0.6 


0.0797783 


1.1 


0.115734 


1.6 


0.184141 


1.99 


0.655754 


0.2 


0.0592869 


0.7 


0.0858022 


1.2 


0.125393 


1.7 


0.209986 


1.995 


0.796158 


0.3 


0.0639583 


0.8 


0.0923055 


1.3 


0.136425 


1.8 


0.248441 


1.999 


1.22703 


0.4 


0.0689056 


0.9 


0.0993773 


1.4 


0.149299 


1.9 


0.321076 


1.99999 


3.98564 



This formula was reported in Ref.|32|for the 5* = 5 antiferro- 
magnetic Heisenberg spin chain. 



B. For finite U 

Next we examine the generic < [/ 7^ 00 case. In this 
case, there appears the conventional 2k p oscillation term in 
addition to the Akp one. We show that the additional loga- 
rithmic correction appears near the phase boundary to the CO 
phase, not only in the Akp oscillation term but in the 2fci? os- 
cillation term. 

Based on the conventional bosonization for electron sys- 
tems, the density operator is given hy^^^ 



p{x) 



1 d(j)p 

TT dx 
. 2C2 



2ci 



sin(2fci?a; 



cos{4:kpx + 20p), 



(8) 



where 0p and 



are the charge and spin phase fields. The ci 
and C2 are nonuniversal numerical quantities satisfying ci ~ 1 
and C2 = in the noninteracting case. In the similar way to 
the U = 00 case, the most general form of the density-density 
correlation function is derived as (see Appendix) 



N{x) 



A. 



-4i 



cos2kpx In ^^^{x/acr) 

2-8i<'„l 1/2 



iAkpx [1 - (a/x) 



[1 + {a/x) 



2-8K, 



i3/2 ' 



(9) 



where a and are the short-distance cutoff parameters for 
the charge and spin sectors, respectively. In the case of 
f/ ^ 00, the coefficients Ai and A2, which are propor- 
tional to Ci and respectively, are to be determined nu- 
merically. The logarithmic correction \n~'^^^{x/acr) in the 
2kF oscillating term appears due to the marginally irrele- 
vant coupling of the spin channel In the noninteracting 
limit U = V = 0, the quantity A2 vanishes and the log- 
arithmic correction ln~^''^(a;/acr) is replaced by a constant, 
and then the correlation function reproduces the trivial result 
N{x) = -l/{-K^x^) + cos2fcFa;/(7r2x2). 

From Eq. (|9]l, we find that an anomalous logarithmic cor- 
rection also appears in the 2kF oscillating term near the 
boundary of the CO phase. On the phase boundary, the corre- 



lation function reads 
1 



N{x) 



47r2x2 



A 



cos Akpx 
' (27r)3/2 



\n'^'{x/c 



for {U,V) ^ (U,,V,) (10) 
where ^1 ^ and A2 ^ I for {Uc, K) ^ (00, 2). 

III. NUMERICAL RESULTS 

For numerical confirmation of the logarithmic corrections, 
we employ the DMRG technique which provides very accu- 
rate data for the ground-state correlation functions of ID cor- 
related electron systemsi^ We consider L/2 electrons on a 
chain with L sites and calculate the equal-time density-density 
correlation function 



N{x) 



(11) 



under the open-end boundary conditions (OBC). Here, the dis- 
tance X is centered at the middle of the system. The appli- 
cation of OBC enables us to obtain the correlation function 
(fTTT i quite accurately for very large finite-size systems up to 
^ 0(1000) sites. However, the real-space DMRG method 
works with finite number of sites, so that we have to pay spe- 
cial attention to the finite-size effects for a precise compari- 
son with the RG results. In the present calculations, the most 
problematic finite-size effect is the Friedel oscillation starting 
from the open edges. To eliminate it, we simply add on-site 
potential energy V/2 on both edge sites. It corresponds to a 
compensation of the "missing correlation" caused by the ab- 
sence of their neighboring site. Hereby the Friedel oscillation 
is fairly suppressed. 

On that basis, the remaining finite-size effects are investi- 
gated. We now choose some parameters in the vicinity of the 
CO phase where the finite-size effect is relatively large due 
to strong charge fluctuations. For these parameters, we cal- 
culate N{x) for several chains with length up to L = 1024 
sites and then obtain the extrapolated values to the thermody- 
namic limit (L 00) using the finite-size-scaling analysis. 
By comparing the extrapolated values and the finite-size data, 
we find that N{x < 200) in the thermodynamic limit can be 
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FIG. 1 : (Color online) (a) DMRG results of the correlation function 
N{x) for {U,V) — (oo, 1.95t). The solid line denotes a fitting 
with Eq. Inset: \N{x) \ plotted on a log-log scale, (b) Estimated 
length scale of the logarithmic correction Xaoaa as a function of V/t 

for U = oo. 
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reproduced with extracted central 200 sites of a chain with 
L = 512 sites within a few percent error The relative error 
in the ground-state energy, |{e(oo) — e(512)}/e(cK))|, is be- 
low 0.1%, where e{L) is the ground-state energy per site for 
a chain with L sites. Consequently, we will study the equal- 
time correlation function N{x) for the central 200 sites of a 
chain with L ^ 512 sites without the finite-size-scaling anal- 
ysis. We keep up to m « 4000 density-matrix eigenstates 
in the DMRG procedure and all the calculated quantities are 
extrapolated to the m ^ oo limit. For your information, in 
this way we obtain (^1,^2) = (0.991,-0.0003) [the exact 
are (^1,^2) = (I7O)] for the coefficients of Eq. Q in the 
non-interacting case U = V = Q which poses a non-trivial 
problem to the DMRG method. 



FIG. 2: (Color online) DMRG results of the correlation function 
N(x) for several sets of {U, V). The solid lines denote fitting curves 
with Eq. l|9ll. The estimated values of Across and Kp are also in- 
cluded. 



are used to carry out the fitting procedure and, however, the fit- 
ting results are almost unchanged for any choice of the length 
from 40 to 300. In the same way, we can also estimate the 
values of Xcioss for the other V/t values. In Fig.[T](b), the 
estimated values of Xcmss are plotted as a function of V/t. 
We find that the logarithmic correction is hardly present at 
V ^ 1.8 and the length scale Xaoss increases rapidly in the 
vicinity of the CO insulating phase. Note that Xcross 00 at 
V/t -> 2. 



A. The f/ -> 00 limit 



B. For finite U 



Let us first consider the correlation function N{x) in the 
U ^ 00 limit. We now attempt to fit the DMRG results of 
N{x) into the analytical form of Eq. (|4]l. Since the exact so- 
lutions of A2 and Kp are available, the quantity a is the only 
fitting parameter in Eq. (|4]i. Figure [T] (a) shows the DMRG 
results of N{x) for ([/, V) = (00, 1.95t). An excellent agree- 
ment of the DMRG data with the fitted line is found. We 
then obtain a — 0.0515, which leads to Xcross = 34.8. It 
means that the logarithmic correction appears at x < 35 for 
V = 1.95t. We note that the central 200 sites out of i = 512 



We next turn to the case of f/ 7^ cxj. In this case, the numer- 
ical results of the correlation function N{x) can be fitted with 
the analytical form of Eq. (|9]l. Differently from the case of 
U = 00, there are five fitting parameters; namely, a, oi„, A\, 
A2, and Kp. Of them, Kp may be obtained very accurately 
with the DMRG method via the derivative of charge structure 
factor at g = 0, 



(12) 



5 



0.4 




FIG. 3: (Color online) Correlation amplitude A2 for [/ = 1 (trian- 
gles), 5 (squares), and 00 (circles). Solid (dashed) lines denote the 
DMRG (exact) results. 
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Thus, we can 



with q = 271 /L and n{q) = ^ 

reduce the fitting parameters from five to four In Fig. |2] we 
show the fitting results of the correlation function N{x) near 
the boundary of the CO phase for U ^ 5t and lOt [the crit- 
ical boundary has been estimated as Vc ~ 3.70t (2.76t) for 
U = 5t (lot) in Ref. HI. We can see that the DMRG data 
is in good agreement with the fitted line for all the parame- 
ter sets. From the obtained results of a;cross, we find that the 
logarithmic correction appears for ([/, V) ~ {5t, 3.5t) and 
(m, 2M). Especially at (U, V) = {5t, 3M), the length scale 
is extremely large Xaoss ~ 1600; it allows us to crossly notice 
that this point is very close to the boundary of the CO phase. 
Meanwhile, the logarithmic correction is hardly present for 
([/, V) = (5<, 3i) and {U, V) = (lOt, 2t). As aresult, we con- 
firm that the logarithmic correction is present also for U ^ 00 
and its length scale grows rapidly near the CO phase bound- 
ary. 

Finally, we discuss the correlation amplitudes, Ai and A2, 
of Eqs. (HI and (|9]l. Figure [3] shows the DMRG results of 
the amplitude A2 as a function of V/t for several values of 
U /t. In the U = 00 limit, we can see an excellent agreement 
between the DMRG and exact results. We also find a very 
sharp increase near V ^ Vc ^ 2t. For U = 5t, the behavior 
of A2 seems to be quite similar to that for U = 00; while, the 
amplitude Ai is rapidly decreased near V ^ Vc = 3.70t, e.g., 
Ai < 10^^ at F > 3. Thus, the 2kp oscillating term would 
be negligible in the vicinity of the CO phase. When U — t, 
the amplitude A2 decreases with increasing V/ 1 in reflecting 
that the Akp fluctuation is not enhanced by V. We note that 
the two amplitudes Ai and A2 are rather small with the same 
order of magnitude for small U and larger V values. 



IV. SUMMARY 

We study the density-density correlation function in the 
TLL state for the ID extended Hubbard model at quarter fill- 
ing. Based on the bosonization and RG techniques, we ob- 
tain the generalized analytical form of the correlation function 



which exhibits anomalous power-law behavior with logarith- 
mic corrections near the phase boundary to the CO insulating 
state. Using the DMRG method, we confirm the appearance 
of the logarithmic corrections not only in the U = 00 limit 
but also for finite U. Moreover, we find that the length scale 
of the corrections grows rapidly near the CO phase boundary. 
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APPENDIX A: DERIVATION OF THE ANALYTICAL FORM 
OF THE CORRELATION FUNCTION 

In this appendix, we derive the generalized form of the cor- 
relation function [Eqs. and (|9|l] based on the RG approach. 
In our derivation, we follow the formalism of the RG method 
developed in Ref. [l^ 

The RG equations for the TLL parameter Kp{l) and the 
1/4-filled umklapp scattering G^H) are given hy^sM 

j^Kpil) = -2G?,(0^p(0, (Ala) 

^Gu(0 = [2 ~ 8A'p(0]Gu(0, (Alb) 

where the initial values are estimated based on the perturbative 
treatment in Ref.|3^ The TLL parameter in the low-energy ef- 
fective theory can be evaluated from the fixed point value of 
Kp{l), i.e., Kp — Kp{oo). The correlation functions for the 
2kp and Akp oscillation parts, defined as C2kF{^' ~ = 
2{sm{2kpx + (f>p{x)) sm{2kpx' + (pp{x'))) / cos 2kp{x — x') 
and C4i;p(x — x') = 2{cos{Akpx + 4)p{x)) sm{2kFx' + 
(f)p{x'))) / cos2kp{x — x'), respectively, are given in the RG 
scheme by2ii2 



C2kF{x) = CXp 

C^kpix) = exp 



ln(x/Qo) 



dlKp{l) 



(A2a) 



\n(x/aQ) 



dl{4Kp{l)-2Guil)) 



(A2b) 



where is the short-distance cutoff. The couplings Kp{l) 
and Gu(0 are determined by solving Eq. dAlb . 

From Eq. (lAlb . we find that the fixed point values are given 
by (A'p(oo), Gu(oo)) ~ on the phase boundary be- 

tween the TLL and CO states. Near this phase boundary, the 
TLL parameter can be expanded as Kp{l) = \ + \Gp{l) and 
we can treat Gp{l) perturbatively. Up to the second order in 
Gp{l) and Gu{l), the RG equations (lAll i are rewritten as 

j^GpH) = -2Gl{l), j^G^{l) = -2Gp(0G,(/). (A3) 

In the case of Gp(0) > |Gu(0)|, the umklapp scattering G^{1) 
flows to zero, i.e., is irrelevant, and Gp{l) has a finite fixed 



6 



point Gp(oo) > 0. Thus the TLL parameter in the low-energy 
limit is given by Kp = \ + \Gp{oo). The explicit solutions 
of Eq. (IA3l l are given by 

Gpil) = -coth[6'/ + tanh"^(6'/2G'p(0))] , (A4a) 

n 

Gu(0 = -cosech[6'/ + tanh"i(6'/2Gp(0))] , (A4b) 

where 9 = 2(Gp — G^)^/^ is a scaling invariant quantity. Near 
the phase boundary, i.e., for small 0, the umklapp scattering 
Gu(0 approaches to zero very slowly as increasing I. By sub- 
stituting Eq. ( IA4I) into Eq. (IA2I) . we obtain the analytical form 
of the correlation functions: 

/ an \ 1/4+8/8 / 1-^28 \ 1/8 

^''^^'^ = It J [ l + {dao/.r ) 



where d is the nonuniversal quantity depending on 
the initial values of RG equations, defined by d = 
cxp[— 0^1 tanh^^(0/2Gp(O))]. In terms of Kp, the param- 
eter 6 is given by = — 2). By defining a = dao we 
can derive Eqs. ^ and (|9|l. 
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